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LONG-TERM  GOAL 

Our  long  term  goal  is  to  develop  multifrequency,  high  frequency  (HE)  radar  techniques  and 
instrumentation  for  measuring  surface  currents,  vertical  shear,  winds,  friction  velocity  and  waves  in 
coastal  regions  and  large  lakes  for  scientific,  civil  and  military  applications.  Our  goal  includes 
deployment  of  HE  radar  systems  for  air-sea  interaction,  coastal  oceanography  and  ship  detection 
research  as  well  as  the  integration  of  HE  radar  measurements  into  coastal  ocean  models.  A  related  goal 
is  to  investigate  and  develop  ship  detection  and  tracking  techniques  for  multifrequency  HE  radar. 

OBJECTIVES 

The  objectives  of  this  project  began  with  the  construction  and  deployment  of  two  and  later  three 
multifrequency  HE  radar  instruments  (called  MCR  for  Multifrequency  Coastal  Radar)  to  Monterey 
Bay,  California.  Eurther  experiment  deployments  were  to  the  Virginia  coast  and  to  Eake  Michigan  for 
fresh  water  experiments  (NSE  sponsorship).  The  data  collected  at  these  sites  is  reduced,  analyzed  and 
interpreted  to  achieve  the  specific  research  objectives  listed  below: 

1.  Improvement  of  radar  performance  by  upgrading  hardware  and  software  and  developing 
improved  transmit  antennas,  signal  processing  and  flexible  use  of  multiple  frequencies 

2.  Improvement  of  HE  radar  estimates  of  surface  currents,  vertical  shear,  winds,  friction  velocity  and 
waves  by  improved  estimation  algorithms  that  use  knowledge  of  air-sea  interaction  physics 

3.  Estimate  surface  wind  speed  and  direction  using  multifrequency  HE  radar  measurements  by 
exploiting  air-sea  interaction  physics  and  HE  radar  measurements  of  vertical  current  shear 

4.  Ocean  science  investigations,  including  assimilation  of  HE  radar  data  into  coastal  ocean  models 
for  circulation,  chemical  and  biological  properties  as  well  as  air-sea  interaction  studies 

5.  Use  of  continuing  HE  radar  observations  on  Monterey  Bay  in  observing  ships  to  assess  the 
usefulness  of  multiple  frequency  HE  radar  in  ship  detection  and  tracking. 
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APPROACH 


This  project  requires  an  expert  team  of  engineers  and  oceanographers  from  many  institutions, 
including  Peter  Hansen  and  Lorelle  M  eadows  (University  of  M  ichigan),  Calvin  Teague  (Stanford 
University),  Dan  Fernandez  (California  State  University,  M  onterey  Bay),  Jeff  Paduan  (Naval 
Postgraduate  School)  and  Kenneth  Laws,  Jessica  Drake  and  Steve  Petersen  (University  of  California  at 
Santa  Cruz).  This  team  designed,  constructed,  upgraded  and  deployed  M  CRs  at  sites  on  the  Atlantic 
and  Pacific  oceans  as  well  as  the  G reat  L akes.  We  continually  seek  to  deploy  M  CR  systems  for  air-sea 
and  ship  observations  in  both  short  campaigns  and  long-term  observations.  A  summary  follows: 

1.  Ocean  science  investigations:  participate  in  investigations  of  air-sea  interaction,  including 
momentum  transfer  and  structure  of  the  air  and  sea  boundary  layers  and  the  assimilation  of  HF 
radar  data  into  the  coastal  ocean  circulation  models 

2.  Improve  radar  performance  by  upgrading  radar  hardware  and  software,  repairing  and 
upgrading  the  radar  antennas,  especially  transmit  antennas;  and  comparing  M  CR  with  other  HF 
systems 

3.  I  mprove  H  F  radar  estimates  of  surface  currents,  current  shear,  winds  and  friction  velocity 

by  developing  improved  estimation  algorithms  that  use  air-sea  interaction  physics,  e.g.,  the 
Charnock  relation,  and  continued  work  on  M  USIC,  rootM  USIC  and  beamforming 

4.  E  stimate  surface  wind  vectors  using  multifrequency  H  F  radar  measurements  by  usi  ng 
statistical  techniques,  air-sea  interaction  theory  and  the  H  F  vertical  current  shear  measurements 

5.  Use  the  continuing  HF  radar  observations  on  Monterey  Bay  to  observeships  (especially 
known  ships  at  known  locations)  at  multiple  HF  frequencies,  noting  changes  in  detection  and 
tracking  performance  with  frequency,  range,  direction,  ship  type,  etc. 

WORK  COMPLETED 

1.  Observational  program:  We  strive  to  keep  the  M  CR  units  in  operation  either  at  M  onterey  Bay, 
their  'home  site',  or  employed  in  experimental  campaigns  elsewhere.  During  the  last  year  two 
M  CRs  operated  on  and  near  M  onterey  Bay.  This  includes  monitoring,  maintenance  and 
upgrades. 

2.  Estimation  of  wind  friction  velocity  (u*)  using  multifrequency  HF  radar:  M  ultifrequency 
HF  radar  observations  of  currents  at  different  depths  are  used  to  find  u*  -  see  M  eadows  (2002). 

3.  Evaluation  of  Root-MUSIC  processing  for  HF  radar:  Weassessed  the  useful  ness  of  this 
specialized  form  of  MUSIC  processing  for  linear  HF  radar  antenna  arrays  -  see  Teague  (2002). 

4.  C  omparison  of  M  C  R  and  SeaSonde  measurements  of  surface  currents  on  M  onterey  Bay 

CA:  We  compared  surface  current  measurement  results  of  MCR  and  SeaSonde  (CODA  R  Ocean 
Systems  instrument)  systems  using  data  from  co-located  sites  -  see  Vesecky  et  al.  (2001  & 

2002). 
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5.  Retrieval  of  vector  winds  using  multifrequency  H  F  radar  measurements:  Partial  least 
squares  estimation  to  estimate  vector  winds;  estimates  were  compared  to  buoy  measurements. 

6.  Communication  of  research  results:  This  year  we  presented  eight  papers  atfour  major 
conferences  and  workshops  on  remote  sensing  and  oceanography,  including  IGARSS  and  the 
Radio  Oceanography  Workshop,  sponsored  by  ONR.  One  refereed  journal  paper  was  published. 

7.  Ph.  D.  student  graduated:  We  are  pleased  that  Lorelle  M  eadows  received  her  Ph.  D.  degree  in 
j  une  2002  from  the  U  niversity  of  M  ichigan  (see  M  eadows,  2002  below). 

RESULTS 

This  year  produced  many  new  results  with  the  analysis  of  past  field  observations  as  well  as  collection 
of  new  data.  Below  are  two  examples  that  are  particularly  important  to  the  field.  One  concerns  the 
retrieval  of  friction  velocity  from  multifrequency  H  F  radar  observations  and  comparison  with  buoy 
measurements  and  the  other  concerns  the  prediction  of  vector  winds  using  H  F  radar  data. 

During  the  M  USE  experiment  of  Autumn,  2000  multifrequency  coastal  radars  (M  CRs)  at  M  oss 
Landing  and  Santa  Cruz  on  California's  M  onterey  Bay  made  FI  F  measurements  that  covered  the 
location  of  the  NPS  flux  buoy  near  the  M  oss  Landing  radar  site.  This  buoy  was  capable  of  measuring 
friction  velocity  in  the  air  as  well  as  other  air-sea  interaction  parameters.  As  her  Ph.  D.  thesis  project, 
Lorelle  M  eadows  used  data  from  these  radars  as  well  as  the  flux  buoy  to  develop  and  test  a  technique 
for  retrieving  friction  velocity  from  M  CR  measurements.  The  approach  to  the  friction  velocity 
measurement  begins  by  assuming  that  the  vertical  current  profile  follows  the  logarithmic  (law  of  the 
wall)  theory  for  a  turbulent  boundary.  Although  there  is  some  controversy  regarding  a  linear  variation 
close  to  the  surface  and  experiments  to  support  both  the  linear  and  the  logarithmic  views,  the 
difference  is  not  critical  to  this  method.  The  fundamental  di rect  measurement  of  the  FI  F  radar  is  the 
Doppler  shift  and  hence  the  phase  velocity  of  the  radially  traveling  Bragg  resonant  waves  on  the  ocean 
surface,  i.e.,  ocean  waves  at  wavelengths  of  ~  6.9, 11.2,  22.1  and  31.3  m.  The  Bragg  resonant  waves 
traveling  radially  with  respect  to  the  radar  are  half  the  radar  wavelength  and  give  rise  to  the  FIF  radar 
echoes.  The  phase  velocity  (Cp)  of  surface  waves  affected  by  near  surface  currents  can  be  estimated  to 
first  order  by  a  weighted  integral  over  depth  z  of  the  vertical  current  profile  U(z),  Cp  =(g/k)  -Fdc 
where  dc  =  2k  jU(z)  e^^^  dz,  where  k  is  the  wavenumber  of  the  Bragg  resonant  ocean  wave  and  g  =  9.8 
m/s^.  The  law  of  the  wall  says  that  U(z)  depends  on  friction  velocity  in  the  water  u*w  and  roughness 
length,  Zq. 

The  Levenberg-M  arquart  nonlinear  estimation  technique  is  used  to  find  values  of  u^w  and  Zo  which 
minimize  the  differences  between  the  observed  and  predicted  values  of  Cp.  Once  u+w  is  estimated,  u*a 
in  the  air  can  be  found  knowing  thatu*  scales  with  density,  u*w  =  (ra/  rw  )°'^  u+a.  Thus,  FIF  estimates 
and  buoy  estimates  of  friction  velocity  can  be  compared.  Note  that  this  method  does  not  rely  on  an 
effective  depth  estimate  for  H  F  radar  currents,  but  uses  the  Doppler  observations  directly. 
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Fig.  1.  Comparison  of  estimates  of  u*a-  The  colored  square  and  diamond  are  estimates  from  the 
NFS  flux  buoy  using  a  bulk  parameter  estimation,  using  air-sea  measurements,  and  a  direct 
calculation  (DC),  found  by  direct  correlation  of  wind  field  fluctuations.  The  multifrequency  coastal 
radar  (MCR)  estimates  are  made  as  described  in  the  text.  Radial  refers  to  along  the  HF  radar  line 
of  sight,  i.e.,  all  the  u*  estimates  are  projected  onto  this  line  of  sight. 


For  wind  speeds  above  ~  6  m/s  the  H  F  radar  estimates  agree  very  well  with  the  buoy  measurements. 

At  lower  wind  speeds  one  expects  poorer  results  from  the  MCR  because  other  factors,  such  as  tides 
and  buoyancy  flows,  are  affecting  the  near  surface  currents  and  are  presently  ignored  in  the  analysis. 

At  any  wind  speed  this  method  produces  values  of  which  tells  one  how  well  the  model  fits  the 
measurements,  and  hence  allows  one  to  put  error  bars  on  the  measurements.  Note  that  aside  from  one 
case  the  buoy  measurements  are  w  ithi  n  the  error  bars  of  the  M  C  R  esti  mates.  We  regard  these  results 
as  very  promising  for  the  use  of  multifrequency  coastal  radar  in  mapping  friction  velocity. 

M  ultifrequency  H  F  radar  measurements  reveal  the  vertical  current  shear  near  the  ocean  surface.  Thus, 
weanticipate  that  wind  speed  and  direction  can  be  estimated  from  MCR  data  since  wind  stress  on 
surface  is  usually  the  dominant  factor  in  determining  current  magnitude  and  the  vertical  shear. 

Previous  research  has  been  successful  in  making  wind  direction  estimates,  but  no  really  successful 
estimates  of  wind  speed.  In  this  study  we  did  not  try  to  employ  a  physical  model  as  above,  but  for  this 
initial  effort  choose  a  linear,  empirical  prediction  method  known  as  partial  least  squares.  This  method 
seeks  a  wind  component  esti  mate  based  on  a  weighted  sum  of  all  the  input  data  components.  The 
weights  are  determined  by  a  least  squares  procedure.  For  the  input  data  we  used  the  surface  currents  at 
all  depths  aswell  asthe  wind  direction  estimated  by  the  relative  strengths  of  the  first  order  Bragg 
peaks  in  echo  Doppler  spectrum.  Weights  derived  by  the  partial  least  squares  algorithm  emphasized 
the  following: 
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1 .  Input  data  fiom  the  highest  radar  frequencies,  sensing  nearest  the  smface,  were  most  important. 


2.  Differences  between  cuixents  estimated  from  radar  data  at  two  different  frequencies  were 
important  -  this  reflects  the  importance  of  shear  in  the  wind  estimates. 

3.  Dir  ection  estimates  fr  om  Bragg  line  ratios  were  important  -  this  ratio  estimates  the  wind 
direction  as  established  by  a  nirmber  of  investigators  in  earlier  work. 


The  results  are  impressive  with  the  standard  enor  of  rneasiuement,  i.e.,  the  standard  deviation  of  the 
difference  between  HF  and  biroy  measirrements,  being  very  near  1  m/s  after  removal  of  a  very  small 
bias,  ~  10  cm/s.  This  level  of  estimation  accruacy  is  soirght  by  meteorologists,  brrt  seldom  attained  in 
remote  sensing  measmements.  Fruther,  the  method  does  not  appear  to  degrade  at  low  wind  speeds. 
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Fig.  2.  Comparison  of  wind  speed  estimates  from  MCR  with  buoy  winds.  The  standard  error  of 
prediction  SEP  is  ~  1  m/s  and  bias  ~  10  cm/s.  Left  (right)  panel  shows  eastward  (westward)  U  (V) 
component  of  the  wind.  Colored  points  are  MCR  estimates,  whereas  the  gray  points  are  the  training 
set  used  to  develop  the  partial  least  squares  estimation  algorithm. 


IMPACT/APPLICATION 

Multifr  equency  HF  radars  have  demonstrated  then  irsefrrlness  hi  measining  smface  ciments  and 
crm  ent  shear  in  the  top  few  meters  of  the  ocean.  No  other  techniqire  makes  sitch  measurements  over  a 
large  area,  at  sirch  low  cost.  We  extended  HF  radar  applications  to  fresh  water  and  to  mapping  of 
fr  iction  velocity  and  wmd  vector.  HF  radar  measmements  are  irsefrrl  for  real  time  applications. 

RELATED  PROJECTS 

1.  We  participated  in  the  Integrated  Coastal  Ocean  Network  (ICON)  project  (National  Ocean 
Partnership  Program),  integrating  in  sitir  and  remote  measmements  over  Monterey  Bay. 
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2.  We  began  participation  in  the  Center  for  I  ntegrated  M  arine  Technology  that  brings  together 
physical,  chemical  and  biological  oceanographers  as  well  as  marine  biologists  to  make  and 
interpret  measurements  relevant  to  the  M  onterey  Bay  M  arine  Sanctuary  (NOAA) 
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